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Abstract 

The extended SSH model and Bogohubov-de Gennes(BdeG) formahsm are 
apphed to investigate the electronic properties and stable lattice configura- 
tions of C36. We focus the problem on the molecule's unusual Dg/j symmetry. 
The electronic part of the Hamiltonian without Coulomb interaction is solved 
analytically. We find that the gap between HOMO and LUMO is small due 
to the long distance hopping between the 2nd and 5th layers. The charge 
densities of HOMO and LUMO are mainly distributed in the two layers, that 
causes a large splitting between the spin triplet and singlet excitons. The 
differences of bond lengths, angles and charge densities among the molecule 
and polar ons are discussed. 
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I. INTRODUCTION 



Recently, a new member of fullerenes, C36, was synthesized by arc-discharge method and 
purified in bulk quantities Up to now, it is the smallest fuUerene ever discovered. The 
C36 molecule is more curved than Ceo because of the adjacent pentagonal rings and the small 
number of carbon atoms [|I|] . This feature suggests stronger electron-phonon interaction and 
possible higher superconducting transition temperatures in the alkali-doped C36 solids than 
those of Ceo 0- The Solid-State Nuclear Magnetic Resonant experiment suggested that 
the most favorable configuration of Cse molecule has Dq^ symmetry This confirmed 
the results of the early ah initio pseudopotential density functional calculations [H, which 
indicated that the D^h structure is one of the most energetically favorable structures among 
several possibilities. 

Looked in a direction perpendicular to the six-fold axis(Fig.l), the molecule is composed 
of six parallel layers. On each layer, six carbon atoms lie at the vertexes of a hexagon. The 
edges of the hexagons of the 1st, 2nd, 5th, 6th layers are parallel, while those of the two 
hexagons in the middle two layers are turned by an angle of 30°. In addition to this special 
structure, the six-fold principal axis is distinct form the ordinary five-fold principal axis in 
the cases of Ceo aiid C70 and little attention has been paid to it. In this paper, we discuss 
both the electronic and structural properties of Cae molecule, emphasizing the interesting 
properties brought by this unusual D^h symmetry. 

We employ a simple and elegant model , extended SSH model, which is successful in 
treating Ceo and C70 [§,11. The BdeG formalism is performed to obtain the stable lattice 
configuration and the corresponding electronic states self-consistently. The electronic part of 
Hamiltonian without Coulomb interaction can be solved analytically by methods of group 
theory. We find that the electron densities of HOMO(-Bi„) and LUM0(-B2g) are zero at 
the two middle layers. Furthermore, the HOMO and LUMO are mainly confined in the 
2nd and 5th layers with a possibility of 90%. The gap between them is considerably small, 
compared to that of Ceo and C70, because the splitting of HOMO and LUMO is due to 
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FIG. 1. The geometrical figure of C36 molecule. 

a long-distance hopping in our case. We find the large sphtting of the triplet and singlet 
excitons due to more localized HOMO an LUMO in comparison with Cgo and C70. As a 
result, the triplet exciton's energy is very small and possible experiments are suggested to 
test this phenomenum. 

The geometrical figures of C36 molecule and polarons are discussed and bond lengths 
and angles are calculated. The charge density on each layer are given. We find that the 
differences of charge densities among polarons and molecule are mainly in the 2nd and 5th 
layers. 

The following sections are arranged as: in Sec.I, the extended SSH model is introduced; 
in SecM, we analytically solve the electronic part of the Hamiltonian without Coulomb 

interaction; in Sec. IV, the whole Hamiltonian is solved sclf-consistcntly; in Sec.V, results 
and discussions are presented; and conclusions are made in the final section. 

II. MODEL 

The Hamiltonian of the extended SSH model to the C36 molecule is written as: 

H — Ho + Hint + Helas 

The first term Hq of Eq.(l) is the hopping term of 7r-electrons. 
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FIG. 2. The geometrical figure of C36 molecule. 

Hq= -^[to- a{lij - lo)]{c\^Cja + h.c.) - Y^ti{c\„Cja + h.c.) 

(iJ) (hj) 

- 12"H4aCja + h.c.) (2) 
(id) 

where the three terms describe hopping between the nearest neighbors, the next nearest 
neighbors and the third neighbors respectively, with the corresponding hopping integrals 
^0,^1 and ^2- The influence of electron- phonon coupling is only included in the first term of 
Hq, since those of the last two terms are much smaller than tg. The third term is important 
here because it changes the accidental degeneracy of HOMO and LUMO states , which will 
be explained in the next section. The larger distance terms are neglected, as they are small 
and do not bring any new effects. 

The second term of Eq.(l) is the screened Coulomb interaction expressed in the Hubbard 
model. 

Hint = Uj^rii^nii + V ^ianja' (3) 

i (ij),<T,<T' 

where U is the strength of the on-site interaction, and V is that between the nearest neigh- 
bors. 

The third term of Eq.(l) is the elastic energy of the lattice. This term is composed of 

three parts, 

Helas = lK,J2ihj - k? + \K^Y.dOi^ + \K^Y.dOi,\ (4) 
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where the first part describes the spring energy of bond-length terms, and the next two terms 
describe the spring energy of the angular terms. = 1 ~ 3) are the elastic constants 

for these different kinds of lattice vibrations. dOi^^dOiQ are bond angle deviations from the 
original angle 108°, 120°. e.g. dOi^ = — 108°. The first summation is taken up over 
all nearest pairs. The second is of all interior angles of pentagons. And the third is of all 
interior angles of hexagons. 

Since there are no enough experiment data to determine the semiempirical parameters, 
we set the values in the scope of fuUerenes such as Ceo and C70 which can produce reasonable 
results. We adjust them to give the gap between HOMO and LUMO and the bond length 
consistent with those of the pseudopotential density functional approach. We take to = 2.5 
eV, a = 5.6 eV/A, Lq = 1.55 A, Ki = 47 eV/ A^, K2 = 8 eV/rad^, = 7 eV/rad^, 
ti = 0.168i(:o- In fact, the particular values, except for the value of t2, would not affect the 
physics too much. Based on the calculation of the third-nearest yr-orbit integral, t2 was 
estimated to be O.lllto- 

III. ANALYTICAL RESULT OF THE ELECTRONIC HAMILTONIAN 
WITHOUT COULOMB INTERACTION 

Fully exploiting the high Dq^ symmetry, we solve Hq algebraically as in the case of Ceo |6|,^. 
For simplicity, we temporarily ignore ti and t2 terms and will discuss their effects in detail later. 
The thirty-six vr orbits form a 36x36 representation of the D^h group. It can be reduced into the 
sum of following irreducible representations: 

{3Aig © Big © 2B2g © SEig © 3E2g} © {2Biu © 3A2u © B2u © SEiu © 3^2«} (5) 

The even-order axis Ce brings properties different from those odd-order axis C5 characterized in 
Ceo and C70, in the way that it has 1-D representations of kind B. 

We reduce this problem by Dq^s subgroup C^h, using the quantum number m(=0,±l,ib2,3) 
which corresponds to the irreducible representation of Ceh, and P{= ±1) to parity. The thirty-six 
vr orbits are recombined as: 
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TABLE I. The analytical solution of Hq without ti and t2 terms. 



m,p 


E(eV) 








0,+(^i9) 


-9.51 


0.581 


0.574 


0.578 




-5.74 


0.791 


-0.232 


-0.566 




5.57 


0.190 


-0.786 


0.589 


0,-{A2u) 


-8.43 


0.816 


0.504 


0.284 




-3.53 


-0.565 


0.589 


0.578 




8.43 


-0.124 


0.632 


-0.765 


±1,+ {Elg) 


-6.28 


0.647 


0.654 


0.340 ± 0.196i 




-1.60 


0.744 


-0.428 


-0.445 ^ 0.257i 




7.66 


0.145 T 0.084 i 


-0.540 ± 0.312 i 


0.763 


±1, - {Elu) 


-8.09 


0.329=f0.190 i 


0.534T0.308 i 


0.690 




-3.25 


0.882 


-0.015 


-0.409^0.236 i 




4.97 


-0.281 


0.787 


-0.476T 0.275 i 


±2,+ {E2g) 


-5.44 


0.096=f0.169 i 


0.289T0.501i 


0.792 




0.269 


0.573 


0.589 


-0.285^0.494 i 




5.40 


0.796 


-0.565 


0.108±0.187 i 


±2, - {E2u) 


-3.07 


0.379 


0.819 


0.216±0.374 i 




3.20 


0.727 


0.026 


-0.343^0.595 i 




6.25 


0.287T0.496 i 


-0.287± 0.497 i 


0.585 


3,+ {B2g) 


-1.13 


0.357 


0.934 





{Big) 


3.08 








1 


{B2g) 


7.73 


0.934 


-0.357 





3, - (i?2n) 


-3.08 








1 


{Blu) 


-1.13 


0.357 


0.934 





{Blu) 


7.73 


0.934 


-0.357 
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\^%) = Y.v'^'{\l,i) + P\7-l,i)}. 



(6) 



where \l,i) represents the iih.{i = 1 ~ 6) carbon atom's tt orbit in the Zth layer. Here I is only 
ranged for 1 to 3, and rj = e^"^!"^ 



The energy eigenstate wavefunction can be expanded with these base vectors. 



rap / J ^mp,l 



mp 



(7) 



=1,3 



Consequently, Hq can also be reduced to 3x3 matrices in the subspace which is spanned by the 
new bases: 



mp 



(8) 



-tb -tc{l + v'"". 

-tc{l + J?"") -tdPrj^"" 

Because of the Dq^ symmetry, there are only four different kinds of bond lengths {La ~ L^) , see 
Fig.2. The corresponding nearest hopping integrals are ti = to — a{Li — Lq) (i = a ~ d) . The energy 
eigenvalues E^^^ are determined by 



+ AX^ + B\ + C = 



where 



rmr 



A = 2ta cos — + {-l)"'pt - d, 



rmr , 



5 = 2cos —[{-l)^ptatd - ti] - {2ti + tt), 
^ , rmr ,^ rmr. 9 , 9 

C = -4C0S ^(1 + cos —)tatl - {-irptlU 



Equation(9) can be solved analytically, so does the eigenvectors. Here we would not list the 
complicated analytical results. Instead, we present numerical results by using parameters sets,ta = 
3.30eV, tb = 2.95eV, tc = 3.24eV, t^ = 3.08eV, which are determined by the bond lengths Lj. The 
coefficients c^^ ; of the energy eigenstate wavefunction and their irreducible representations 
of Deh that they belong to are showed in Table 1. 

In the C36 molecule, because the hexagons in two middle layers turn an angle of 30° 
to the four hexagons at two ends, the wavefunctions of and Big energy eigenstates 



are entirely composed of the atom orbits on the two middle layers, while those of Biu 
or energy eigenstates only have amplitude on the other four layers. There are six 
energy eigenstates of the B kind representations. The 13th(i?2u), 22nd(i?ip), the accidentally 
degenerate 35th(Siu) and 36th(B2s) levels lie either far below or above the fermi level, so 
they are of no interest. There are also another pair of accidentally degenerate 18th(i?i,j) 
and 19th(i?25) levels which are half filled. We also notice in Table 1. that the 18th(i?iu) and 
19th(S2g) electrons are distributed in the 2nd and 5th layers with 90% possibility and 10% 
in the 1st and 6th layers. 

When ti(16.8%to) and t2(ll-l%^o) are considered, H^p have a different form. The long 
distance hopping(i2) coupling the 2nd and 5th layers results in and -B2s's splitting. This 
is the reason that we introduce much longer distance hopping term than usual. Because 
the charge densities of B^^ and B2g are mainly distributed in the 2nd and 5th layers, the 
gap between HOMO(Bi„) and LUM0(S2g) is approximately 2t2. When the electron-electron 
interaction is taken into account, the electronic Hamiltonian has to be solved in the HF mean 
field theory with self-consistent method. The energy levels of ground state configuration are 
shown in Fig. 3. We can see that in Fig. 3 and Table 1, the relative positions of the energy 
levels only change shghtly, except the splitting of Bi^ and B2g. In Fig. 2, the separation 
between the 16th, 17th levels(£'ip) and HOMO is about 0.2eV and that between the 20th, 
21st levels (^2^) and LUMO is 1.3eV. 

IV. NUMERICAL RESULT OF THE MOLECULE 

Under the adiabatic approximation, we apply the BdeG formalism to the Hamiltonian 
of the molecule to obtain the stable lattice configuration and the corresponding electronic 
structure. Molecular dynamical procedure is used to gradually approach the minimum of 
the potential surface from an initial lattice configuration. The following two equations are 
used, 

F.. = -^ (<^ = l-3) (9) 



= ^ (10) 

where is the effective lattice potential that includes the elastic energy of the lattice and 
the electron-lattice interaction energy, both of which depends on the lattice coordinates. F^^ 
is the a component of the effective force acted on the ith atom under this potential, xia is 
the (7 coordinate of the ith atom. 

In each step of the dynamical procedure, we solve the electrical part of the Hamiltonian 
under the lattice configuration given by the last step. The HF mean field theory is performed 
to decouple the electron-electron interaction. Then the electronic states and wavefunctions 
are obtained self-consistently. 

i 

{i,j),a,a' 

The ground state and the lowest triplet and singlet excitons are investigated. There are four 
lowest exciton configurations A, B, C and D, see Fig. 4. Due to Coulomb interaction, C 
and D are mixed to give singlet and triplet. 



-^E{CT«V';iO-)V'^T«V'ai(i)+CTO>;i«^^T(j)V'aiW} (12) 

where a and (3 denote the 18th and 19th energy level respectively. 



V. RESULTS AND DISCUSSION 

The configurations and energy of ground and low excited states are given below, 

Ef = -64.25 eV |<I>i) (spin singlet ground state), (13a) 

Ef = -64.00 eV 1*^^^) = \B),-^{\C) + \D)) (spin triplets), (13b) 

v2 

= -63.37eV, |$3) = i(|C) - l-D)) (spin singlet). (13c) 
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1 


2 


3 


4 


5 


6 


p2- 
^36 

^36 
C36 

r+ 

*-^36 
p2+ 


1.053 
1.042 
1.031 
1.021 
1.013 


1.131 
1.057 
0.981 
0.902 
0.823 


0.991 
0.985 
0.988 
0.992 
0.997 


0.991 
0.985 
0.988 
0.992 
0.997 


1.131 
1.057 
0.981 
0.902 
0.823 


1.053 
1.042 
1.031 
1.021 
1.013 



TABLE II. The ground state charge density per site on each layer of ions and molecule. 

We can see that the ground state is a spin singlet rather than triplet after the long 
distance hopping ^2 and electron-electron interaction are taken into account, which is an 
improvement of the prediction of simple Hiickel theory, and there is no need to add the 
hybridizing of a and vr bonds to investigate the qualitative physics of C36 as discussed in 
Ref. [g. 

The lowest excited states \l/2~^ are spin triplet excitons, which are only about 0.25eV 
above the ground states. The energy of singlet exciton \l/3 is very high, and the splitting 
between the triplet and singlet excitons is about 0.63eV, which is much larger than 0.2eV 
in the cases of Ceo and C70. This is because the electrons of HOMO and LUMO are more 
localized, 90% in the 2nd and 5th layers and 10% in the 1st and 6th layers. Compared to 
the C36 case, the electron densities of HOMO and LUMO in Ceo and C70 are distributed 
more uniformly over all sites, so the splitting of the singlet and triplet exciton due to the 
Coulomb interaction is relatively small. The fact that the triplet exciton has low energy 
can be verified through experiment observation: When illuminated by external light, the 
Cae can be excited to the state of singlet exciton and then may transit to triplet exciton 
through nonradiative decay. The transition rate to ground state is slow because of the 
different spin configuration and small energy splitting. We predict that the triplet excitons 
are metastable states, and the phosphorescence phenomena is possible to be observed. The 
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Li 


L2 


L3 


L4 


SSH(C2-) 


1.414 


1.453 


1.418 


1.447 


SSH (C3-6) 


1.410 


1.462 


1.418 


1.447 


SSH(C36) 


1.407 


1.470 


1.418 


1.447 


SSH(C3+) 


1.404 


1.477 


1.417 


1.447 


SSH(C2+) 


1.402 


1.484 


1.417 


1.448 


LDA(C36) 


1.41 


1.48 


1.43 


1.43 




TABLE III. The bond length(A) of 


ions and molecule. 




a 




7 


6 e 




C36 120.00° 


107.64° 


108.20° 


108.29° 119.42° 


119.24° 



TABLE IV. The bond angle within SSH model. 

triplet exciton is paramagnetic while the ground state is a diamagnetic singlet. So the 
magnetic susceptibility is increased upon illumination of the external light. Electron Spin 
Resonance(ESR) experiment can be performed to detect the triplet exciton states. When 
ESR is performed to a solution sample, the usual unimportant magnetic dipole interaction 
between electrons becomes crucial because its anisotropy can smear the resonance peak. 
However, this interaction is decayed rapidly with distance as R~^. The electrons in HOMO 
and LUMO are mainly uniformly distributed on 12 carbon atoms in two layers, the possibility 
of their short distance is considerably small compared to other small organic molecule, such 
as naphthalene. So the resonance peak is possible to be observed. The phone absorption 
and luminescence spectra experiments can be perform to test the lowest singlet exciton's 
energy. 
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The shape of the C36 molecule is an eUipsoid with high aspect ratios. Our calculation 
shows that the distance between the 1st and 6th layers is 5.2A°, and that between the 
parallel vertical hexagon planes is 4:.2A°. The charge density, bond lengths and bond angles 
of C36 molecule and ions are calculated and shown in Table 2, 3 and 4. We note that because 
of the DQfi symmetry of the molecule, the 1st and 6th layers are regular hexagons, with each 
internal angle 120° and the shortest bonds (Li) 1A1A°. The six vertical hexagons around 
the equator are shghtly deformed due to the drag of the 1st and 6th layers from each side. 
L4 is 2.8% longer than Li and 2% than L3. The L2 bonds are relatively weaker than those 
in hexagons, because they couple the hexagons in the two ends with the hexagons around 
the equator. L2 is 4.4% longer than Li. This result is consistent with that of LDA except 
that the L4 is as long as L3 in LDA, and we think that our result is more reasonable. 

The bond lengths of L3 and L4 of ions differ from those of molecule very slightly. But the 
negative polarons' Li is longer and positive polaron' is shorter than the neutral molecule's, 
while the L2 changes in the opposite way. 

The reason is that the wavefunctions of the electron in Bi^ and has opposite sign 
between two nearest sites connected by Li, same sign by L2, and have zero amplitude on 
the 3rd, 4th layers. According to the contribution f these two states form hopping term, the 
sites connected by Li are repulsive, but those connected by L2 are attractive. The angles in 
the pentagons and hexagons around the equator deviate from 108° and 120° shghtly, which 
justifies the validity of H^ias- The bond angles in polarons differ from those of molecule 
very slightly, so we omit them in Table 4. Form C3g' polaron to C^q , the adding electrons 
distribute mostly in the 2nd and 5th layers. This is because of the particular charge density 
distribution of the 18th and 19t levels. 

VI. CONCLUSION 

In this paper, we have carefully studied the effect of symmetry on Cse's elec- 
tronic properties under the extended SSH model. A small gap between H0M0(B2u) and 
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LUMO(Sip) is obtained due to long distance hopping. The large splitting of the spin triplet 
and singlet lowest excitons, the differences of bond lengthes and electron density between 
molecule and polarons are discussed as results brought by the more localized HOMO and 
LUMO and their special symmetries. Possilbe experiments are suggested. 
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